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We present a data structure that supports three-dimensional range reporting queries in 
0(loglog[7 + (log log n) 3 + k) time and uses 0(nlog 1+e n) space, where U is the size of the 
universe, k is the number of points in the answer, and e is an arbitrary constant. This re- 
sult improves over the data structure of Alstrup, Brodal, and Rauhe (FOCS 2000) that uses 
0(n log 1+e n) space and supports queries in (9(logn + k) time, the data structure of Nekrich 
(SoCG'07) that uses 0(n log 3 n) space and supports queries in 0(loglog[/ + (log log n) 2 + k) 
Q \ time, and the data structure of Afshani (ESA'08) that uses 0(n log 3 n) space and also supports 

queries in 0(log log U + (log log n) 2 + k) time but relies on randomization during the preprocess- 
ing stage. Our result allows us to significantly reduce the space usage of the fastest previously 
known static and incremental d-dimensional data structures, d > 3, at a cost of increasing the 
\ query time by a negligible O(loglogn) factor. 
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The range reporting problem is to store a set of d-dimensional points P in a data structure, so that 
for a query rectangle Q all points in Q n P can be reported. In this paper we significantly improve 
the space usage and pre-processing time of the fastest previously known static and semi-dynamic 
data structures for orthogonal range reporting with only a negligible increase in the query time. 

The range reporting is extensively studied at least since 1970s; the history of this problem is rich 
with different trade-offs between query time and space usage. Static range reporting queries can 
be answered in 0(log d n + k) time and 0(n log d_1 n) space using range trees [4] known since 1980; 
here and further n denotes the number of points in P and k denotes the number of points from P 
in the query rectangle. The query time can be reduced to 0(log d_1 n + k) time by applying the 
fractional cascading technique of Chazelle and Guibas [8] designed in 1985. The space usage was 
further improved by Chazelle [6j. In 90s, Subramanian and Ramaswamy |13| and Bozanis, Kitsios, 
Makris, and Tsakalidis [5] showed that d-dimensional queries can be answered in 0(log d_2 n + k) 
timd3 at a cost of higher space usage: their data structures use 0(n log^ 1 n) and 0(n log d n) space 
respectively. Alstrup, Brodal, and Rauhe [2] designed a data structure that answers queries in 
C^log^ -2 n + k) time and uses 0{n\og d - 2+£ n) space for an arbitrary constant e > 0. Nekrich |12j 
reduced the query time by O(logn) factor and presented a data structure that answers queries 
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define 0(/(n)) = 0(/(n) log c (/(n))) for a constant c. 
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Table 1: Data structures in d > 3 dimensions; f indicates that a data structure is randomized. We 
define log*(n) = min{ t | log^ n < 1 } and log** n = min{ t | log*^ n < 1 } where log*^ n denotes 
computing log* t times. 

in 0{\.og d ~ 3 n/(loglogn) d_5 + k) time for d > 3. Unfortunately, the data structure of [12] uses 
0(n log d+1+£ n) space. Recently, Afshani [1] reduced the space usage to 0(nlog d+£ n); however his 
data structure uses randomization (during the preprocessing stage). In this paper we present a 
data structure that matches the space efficiency of [2] at a cost of increasing the query time by a 
negligible O(loglogn) factor: our data structure supports queries in 0(log d ~ 3 n/(loglogn) d ~ 6 + k) 
time and uses 0(n log d ~ 2+e n) space for d > 3. See Table Q] for a more precise comparison of 
different results. 

Our result for d-dimensional range reporting is obtained as a corollary of a three-dimensional 
data structure that supports queries in 0(loglogC7 + (log log n) 3 + k) time and uses 0(ralog 1+£ n) 
space, where U is the size of the universe, i.e. all point coordinates are positive integers bounded 
by U. Our three-dimensional data structure is to be compared with the data structure of [2] that 
also uses 0{n\og l+e n) space but answers queries in 0(log n + k) time and the data structure of [12] 
that answers queries in 0(log log U + (log log n) 3 + k) time but needs 0(n log 4 n) space. See Table[2] 
for a more extensive comparison with previous results. A corollary of our result is an efficient semi- 
dynamic data structure that supports three-dimensional queries in 0(logn + k) time and insertions 
in 0(log 5 n) time. Thus we improve the space usage and the update time of fastest previously 
known semi-dynamic data structure [12] that supports insertions in 0(log 8 n) time. 

If we are ready to pay penalties for each point in the answer, the space usage can be further 
reduced: we describe a data structure that uses 0{n\og d ~ 2 n(loglogn) 3 ) space and answers queries 
in 0(log d_3 n (log log n) 3 + k log log n) time. We can also use this data structure to answer emptiness 
queries (to determine whether query rectangle Q contains points from P) and one-reporting queries 
(i.e. to report an arbitrary point from P(~)Q ifPnQ ^ 0). This is an O(logn) factor improvement 
in query time over the data structure of Alstrup et. al. [2j. Other similar data structures are either 
slower or require higher penalties for each point in the answer. 

Throughout this paper, e denotes an arbitrarily small constant, and U denotes the size of the 
universe. If each point in the answer can be output in constant time, we will sometimes say that the 
query time is 0(f(n)) (instead of 0(f(n) + k)). We let [a, b] denote the set of integers {i\a < i < b}; 
The intervals [a, b) and (a, b] denote the same set as [a, b] but without a (resp. without b). We denote 
by [b] the set [1,6]. 

In section [3] we describe a space efficient data structure for three-dimensional range reporting 
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Table 2: Three-dimensional data structures; f indicates that a data structure is randomized. 

on a [U] x [U] x [{/] grid, i.e. in the case when all point coordinates belong to [U]. In section [5] we 
describe a variant of our data structure that uses less space but needs O (log log n) time to output 
each point in the answer. All results of this paper are valid in the word RAM computation model. 

2 Preliminaries 

We use the same notation as in |14j to denote the special cases of three-dimensional range reporting 
queries: a product of three half-open intervals will be called a (l,l,l)-sided query; a product 
of a closed interval and two half-open intervals will be called a (2,l,l)-sided query; a product 
of two closed intervals and one half-open interval (resp. three closed intervals) will be called a 
(2,2,l)-sided (resp. (2,2,2)-sided) query. Clearly (l,l,l)-sided queries are equivalent to dominance 
reporting queries, and (2,2,2)-sided query is the general three-dimensional query. The following 
transformation is described in e.g. [14] and [13] . 

Lemma 1 Let 1 < ctj < b{ < 2 for i = 1,2,3. A data structure that answers (01,02,03) queries in 
0(q(n)) time, uses 0(s(n)) space, and can be constructed in 0(c(n)) time can be transformed into 
a data structure that answers (01,62,63) queries in 0(q(n)) time, uses 0(s(n) log* n) space and can 
be constructed in 0(c(n) log* n) time for t={b\ — a\) + (62 — 02) + (63 — 03). 

We say that a set P is on a grid of size n if all coordinates of all points in P belong to an 
interval [n]. We will need the following folklore result: 

Lemma 2 There exists a 0{n l+e ) space data structure that supports range reporting queries on a 
d-dimensional grid of size n for any constant d in 0(k) time. 

Proof: One dimensional range reporting queries on the [n] x [n] x [n] grid can be answered in 
0(k) time using a trie with node degree n £ . Using range trees [1] with node degree p we can 
transform a (i-dimensional 0(s(n)) space data structure into a {d + l)-dimensional data structure 
that uses 0{s{n)h{n) ■ p) space and answers range reporting queries in 0(q(n)h(n)) time, where 
h(n) = log nj log p is the height of the range tree. Since p = n e , h(n) = 0(1). Hence, the query 
time does not depend on dimension and the space usage increases by a factor 0(n £ ) with each 
dimension. □ 

We use Lemma [2] to obtain a data structure that supports queries that are a product of a (d — 1)- 
dimensional query on a universe of size n 1_e and a half-open interval. We will show in the next 
Lemma that such queries can be answered in 0(n) space and O(l) time. 
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Lemma 3 There exists a 0(n) space data structure that supports range reporting queries of the 
form Q' x [—00, x) in 0(k) time, where Q 1 is a (d — \) -dimensional query on [Z7i] x [U2] X . . . x [J7<2-i] 
andUi-U 2 - ■■■■U d - 1 = 0{n l - £ ). 

Proof: There are 0(n 1_e ) possible projections of points onto the first d — 1 coordinates. Let 
min(pi, . . . ,Pd-i) denote the point with minimal <i-th coordinate among all points whose first 
d — 1 coordinates equal to pi,P2, ■ ■ ■ ,Pd-l- We store points min(pi, . . . ,Pd-i) for all p\ 6 [f/i],P2 G 
\Ui\i- ■ - ,Pd-i ^ [^d-i] i n a data structure M. Since M contains 0(n l ~ £ ) points, we can use Lemma[2] 
and implement M in 0(n) space. For all possible p\ G [f7i],£>2 £ [U2] >■ • Vd-i G [Z/<j— 1] we also store 
a list L(pi, . . . ,p d ~i) of points whose first d— 1 coordinates arepi, . . . ,Pd-i) points in L(pi, . . . ,Pd-i) 
are sorted by their <i-th coordinates. Given a query Q = Q' x [—00,2;), we first answer Q using 
the data structure M. Since M contains Ofa 1 " 6 ) points, we can find all points in M n Q in 
0(\M n Q\) time. Then, for every point p = (pi, . . . ,pd~i,p d ) found with help of M, we traverse 
the corresponding list L(p\, . . . ,Pd-i) and report all points in this list whose last coordinate does 
not exceed x. □ 

In several places of our proofs we will use the reduction to rank space technique [TOl E]. 
This technique allows us to replace coordinates of a point by its rank. Let P x , P y , and P z 
be the sets of x, y-, and z-coordinates of points from P. For a point p = (p x ,Py,Pz), let 
p' = (rank(pa;, P x ), XBXxk(p y , P y ), rank(p z , P z )), where rank(e,S') is defined as the number of ele- 
ments in S that are smaller than or equal to e. A point p belongs to an interval [a, b] x [c, d] x [e, /] 
if and only if a point p' belongs to an interval [a',b'] x [c',d'] x [e',f] where a' = succ(a,P x ), 
b' = pred(6, P x ), d = succ(c, P y ), d' = pred(d, P y ), e' = succ(e, P z ), f = pred(/, P z ), and succ(e, S) 
(pred(e, S)) denotes the smallest (largest) element in S that is greater (smaller) than or equal to 
e. Reduction to rank space can be used to reduce range reporting queries to range reporting on 
the [n] x [n] x [n] grid: Suppose we can find pred(e, s) and succ(e, S) for any e, where S is P x , P y , 
or P z , in time f(n). Suppose that range reporting queries on [n] x [n] x [n] grid can be answered 
in time 0(g(n) + k). Then we can answer range reporting queries in 0(f(n) + g(n) + k) time. 
Following [2j, we can also use the reduction to rank space technique to reduce the space usage: if 
a data structure contains m elements, reduction to rank space allows us to store each element in 
O(logm) bits. 

3 Space Efficient Three-Dimensional Data Structure 

In this section we describe a data structure that supports three-dimensional range reporting queries 
in 0((loglogn) 3 + log log U + k) time where U is the universe size and uses 0(nlog 1+e n) space. 
Our data structure combines the recursive divide-and-conquer approach introduced in [2], the 
result of Lemma El and the transformation of {a\, 0,2, a3)-queries into (61, &2> &3)-queries described 
in Lemma [TJ We start with a description of a space efficient modification of the data structure for 
(l,l,l)-sided queries on the [n] x [n] x [n] grid. Then, we obtain data structures for (2, 1, l)-sided 
and (2, 2, l)-sided queries on the [n] x [n] x [n] grid using the recursive divide-and-conquer and 
Lemma [3l Finally, we obtain the data structure that supports arbitrary orthogonal queries on the 
[n] x [n] x [n] grid using Lemma [H Reduction to rank space technique described in section [2] allows 
us to transform a data structure on the [n] x [n] x [n] grid into a data structure on the [U] x[U]x [U] 
grid , so that the query time increases by an additive term O (log log U) and the space usage is not 
increased. 
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Lemma 4 ffflj Given a set of three-dimensional points P and a parameter t, we can construct in 
0(nlog 3 n) time a 0(n) space data structure T that supports the following queries on a grid of size 
n: 

(i) for a given query point q, T determines in 0((log log n) 2 ) time whether q is dominated by at 
most t points of P 

(ii) if q is dominated by at most t points from P, T outputs in 0(t + (log log n) 2 ) time a list L of 
0(t) points such that L contains all points of P that dominate q. 

As described in |12| . Lemma H] allows us to answer (l,l,l)-sided queries in 0((log log n) 2 ) time 
and 0{n log n) space. We can reduce the space usage to O(nloglogn) using an idea that is also 
used in p]. 

Lemma 5 There exists a data structure that answers (l,l,l)-sided queries on [n] x [n] x [n] grid 
in 0((log log n) 2 + k) time, uses 0(n log log n) space, and can be constructed in 0(n log 3 n log log n) 
time. 

Proof: For each parameter t = 2 2i , i = i m m,i m - m +i, ■ ■ ■ ,loglogn/2, z m j n = 21ogloglogn, we 
construct a data structure Tj of Lemma [H Given a query point q, we examine data structures Tj, 
i = ^mini *min+i ; • • • > log logn/2 until q is dominated by at most 2 2i points of P or the last data 
structure Tj is examined. Thus we identify the index I, such that q is dominated by more than 2 2 ' 
and less than 2 2 ' +2 points or determine that q is dominated by at least log re points. If I = i m i n , 
then q is dominated by 0((loglogn) 2 ) points. We can generate in 0((loglogn) 2 ) time a list L of 
0((loglogn) 2 ) points that contains all points dominating q. Then, we examine all points in L and 
output all points that dominate q in 0((log log re) 2 ) time. If log log re/2 > I > z m i n , we can examine 
data structures Tj m , n , Tj min+ i,. . ., 7] in 0((l — i m i n )(loglogre) 2 ) time. Then, we generate the list L 
that contains all points that dominate q in 0(2 2 ') time. We can process L and output all k points 
that dominate q in 0(2 2 ') time. Since k > 2 2/ ~ 2 , k = VL(2 2[ ) and k = - i min ) ■ (loglogre) 2 ). 
Hence, the query is answered in 0(k) time. If / = loglogn/2, then q is dominated by r2(logn) 
points, in this case we can use a linear space data structure with O(logre) query time, e.g. the 
data structure of Chazelle and Edelsbrunner [7], to answer the query in 0(logn + k) = 0(k) time. 

Since each data structure Tj uses linear space, the space usage of the described data structure 
is O (re log log n). □ 

Lemma 6 There exists a data structure that answers (2,l,l)-sided queries on [re] x [re] x [re] grid 
in 0((loglogn) 3 + k) time, uses 0(nlog £ n) space, and can be constructed in 0(n log 3 n log log n) 
time. 

Proof: We divide the grid into x-slices = [xj-i, Xj\ x [n] x [n] and y-slices Yj = [n] x [yj-i, yj] x [n], 
so that each x-slice contains n 1 / 2+7 points and each y-slice contains n 1//2+7 points; the value of a 
constant 7 will be specified below. The cell Cij is the intersection of the i-th x-slice and the j-th y- 
slice, = XiPiYj. The data structure D t contains a point (i, j, z) for each point (x, y, z) G PnCij. 
Since the first two coordinates of points in Dt are bounded by n 1//2_7 , Dt uses O(n) space and 
supports (2,l,l)-sided queries in constant time by Lemma[3j For each x-slice Xi there are two data 
structures that support two types of (l,l,l)-sided queries, open in +x and in — x directions. For 
each y-slice Yj, there is a data structure that supports (1, 1, l)-sided queries open in +y direction. 
For each y-slice Yj and for each x-slice Xi there are recursively defined data structures. Recursive 
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Figure 1: Example of a (2, 1, l)-sided query projected onto the xy-plane. i\ = 1, %2 = 5, j\ = 4 and 
ao = xi, b = x 4 , c = V3- 



subdivision stops when the number of elements in a data structure is smaller than a predefined 

constant. Hence, the number of recursion levels is flog log n for v = log 2 2. 

1+27 

Essentially we apply the idea of [2] to three-dimensional (2, 1, l)-sided queries. If a query spans 
more than one x-slab and more than one y-slab, then it can be answered by answering two (1, 1, 1)- 
sided queries, one special (2, 1, l)-sided query that can be processed using the technique of Lemma[3j 
and one (2, 1, l)-sided query to a data structure with n 1 / 2 " 4 " 7 points. If a query is contained in a 
slab, then it can be answered by a data structure that contains n 1 / 2+7 points. We will show below 
that the query time is 0((log log n) 3 ). Each point is stored in 0(2 l ) data structures on recursion 
level i, but space usage can be reduced because the number of points in data structures quickly 
decreases with the recursion level. We will show below that every point in a data structure on 
recursion level i can be stored with approximately (logn/2*) log £ n bits for an arbitrarily small e' . 

Query Time. Given a query Q = [a, b] x (—00, c] x (—00, d] we identify the indices i±, 12, and 
ji such that projections of all cells C{j, i\ < i < 12, j < ji, are entirely contained in [a, b] x (— 00, c]. 
Let ao = ^0 = ^i 2 -l> an d Co = Uj 1 -i- The query Q can be represented as Q = Q1UQ2UQ3UQ4, 
where Q\ = [a ,b ] x (-00, c ] x (-00, d], Q 2 = [a, ao) x (-00, c] x (-00, d], Q 3 = (b ,b\ x (—00, c] x 
(— 00, d], and Q4 = [ao, bo] x (co, c] x (—00, d]. See Fig. [T]for an example. Query Q\ can be answered 
using D t . Queries Q2 and Q3 can be represented as Q2 = ([— oo,ao) x (— 00, c] x (— 00, d]) n X; Ll 
and Q3 = ((—00, b] x (—00, c] x (—00, d]) n Xi 2 ; hence, Q2 and Q3 are equivalent to (1, 1, l)-sided 
queries on x-slices and Xi 2 . The query Q4 can be answered by a recursively defined data 
structure for the y-slice Yj 1 because Q4 = ([ao,bo] x ( — 00, c] x (— 00, d]) n Yj x . If i\ = 12 and the 
query Q is contained in one x-slice, then Q is processed by a recursively defined data structure for 
the corresponding x-slice. Thus a query is reduced to one special case that can be processed in 
constant time, two (1, 1, l)-sided queries, and one (2,l,l)-sided query answered by a data structure 
that contains n 1//2+7 elements. 

Queries Q2 and Q3 can be answered in 0((log log n) 2 ) time, the query Q\ can be answered in 
constant time. The query Q4 is answered by a recursively defined data structure that contains 
0(n l / 2+ i) elements. If i\ = 12 or j\ = 1, i.e. if Q is entirely contained in one x-slice or one y-slice, 
then the query is answered by a data structure for the corresponding slice that contains 0{n l / 2+ ^) 
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elements. Hence, the query time q(n) = 0((loglogn) 2 ) + g(n 1//2+7 ) and q(n) = 0((loglogn) 3 ). 

Space Usage. The data structure consists of O(loglogn) recursion levels. The total number 
of points in all data structures on the i-th recursion level is 2 l n. Hence all data structures on the 
i-th. recursion level require 0(2*nlogn) bits of space. The space usage can be reduced by applying 
the reduction to rank space technique [lOl E] . As explained in section [2j reduction to rank space 
allows us to replace point coordinates by their ranks. Hence, if we use this technique with a data 
structure that contains m elements, each point can be specified with O(logm) bits. Thus, we can 
reduce the space usage by replacing point coordinates by their ranks on certain recursion levels. 

We apply reduction to rank space on every 51oglogn-th recursion level for 5 = e/3. Let V be 

an arbitrary data structure on recursion level r = s^loglogn — 1 for 1 < s < (1/5) log 2 2. Let W 

1+27 

be the set of points that belong to an x-slice or a y-slice of V. We store a dictionary that enables 
us to find for each point p = (p x ,Py,Pz) from W a point p' = (p' x ,p' y ,p' z ) where p' x = rank(p x , W x ), 
p' y = rank(p y , W y ), p' z = rank(p 2 , W z ), and W x ,W y , and W z are the sets of x-, y-, and ^-coordinates 
of all points in W. Let W be the set of all points p' . Conversely there is also a dictionary that 
enables us to find for a point p' G W the corresponding p £ W. The data structure that answers 
queries on W stores points in the rank space of W. In general, all data structures on recursion 
levels r, r + 1, . . . , r + <51oglogn — 1 obtained by subdivision of W store points in rank space of 
W. That is, point coordinates in all those data structures are integers bounded by \W\. If such 
a data structure R is used to answer a query Q, then for each point pr 6 RC\ Q, we must find 
the corresponding point p € P. Since range reduction was applied 0(1) time, we can find for any 
Pr £ R the corresponding p £ P in 0(1) time. 

Each data structure on level r = s<51oglogn for < s < (1/S)v and v = 1^(2/(1+27)) con " 
tains 0(n l ) elements for / = (1/2 + j) r . Hence an arbitrary element of a data structure on level 
r can be specified with I ■ logra bits. The total number of elements in all data structures on 
the r-th level is n2 r . Hence all elements in all data structures on the r-th recursion level need 
0(n2 r ((±4^) r )lograloglogn) bits. 

We choose 7 so that (1 + 2 7 ) < 2 5 / 2 . Then v = 1-t } > 1^75 an d (1 + 27) < 2 s / 2 < 

2 «5-5 2 /2 < 2 &/v = 2 £/3f_ Since r < „io g logn, (1 + 2 7 ) r < 2( £ / 3 ) lo s lo s n < log e/3 n. Therefore 
all data structures on level r use log £//3 n • 0(n log n log log n) = 0(n log 1+2£ / 3 n) bits of space or 
0{n log 2£//3 n) words of log n bits. The number of elements in all data structures on levels r + 1, r + 
2, . . . increases by a factor two in each level. Hence, the total space (measured in words) needed 
for all data structures on all levels q, r < q < r + 5 log log n, is 

(E/=i 1 2 / )0(nlog 2£/3 n) = 

O ( n2 5iogiogn nlog 2 £ /3 n ) = o{n\og e n) because 5 < e/3 and 2 51o s lo s n < log e/3 n. Thus all data 
structures in a group of 5 log log n consecutive recursion levels use 0(n log e n) words of space. Since 
there are (1/5) v = 0(1) such groups of levels, the total space usage is 0(nlog e n). 

Construction Time. The data structure on level (the topmost recursion level) can be 
constructed in 0(n log 3 n log log n) time. The total number of elements in all data structures on 
level s is 2 s nloglogn. But each data structure on the r-th recursion level contains at most n r = n l 
elements and can be constructed in 0(l 3 ■ n r log 3 n log log n) time where I = (1 + 2 r y) r /2 r . Hence, 
all data structure on the r-th recursion level can be constructed in 0((2 r / 3 )nlog 3 n log log n) = 
0(((1 + 27) 3r /2 2r )n log 3 nloglogn) time. We can assume that e < 1. Since we chose 7 so that 
(1 + 27) < 2 £ / 6 , (1 + 2 7 ) 3 < 2; hence, (1 + 2 7 ) 3 72 2r < l/2 r . Then, all data structure on the r-th 
recursion level can be constructed in 0((l/2 r )n log 3 nloglogn) time. Summing up by all r, we see 
that all recursive data structures can be constructed in 0(n log 3 nloglogn) time. □ 
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Figure 2: Example of a (2,2,l)-sided query projected onto the xy-pl&ne. i\ = 2, 12 = 7, j\ = 2, and 
= 5. 



Lemma 7 There exists a data structure that answers (2, 2,1) -sided queries on [n] x [n] x [n] grid 
in 0((loglogn) 3 + k) time, uses 0{n\og £ n) space, and can be constructed in 0(n log 3 n log log n) 
time. 

Proof : The proof technique is the same as in Lemma [6l The grid is divided into x-slices Xi = 
, Xi]xnxn and y-slices Yj = nx[yj_i,yj]xn in the same way as in the proof of LemmaEl Each 
x-slice Xi supports (2, 1, l)-sided queries open in +x and — x direction; each y-slice Yj supports 
(2, 1, l)-sided queries open in +y and — y direction. All points are also stored in a data structure 
D t that contains a point (i,j,z) for each point (x,y,z) £ Pfl Cy. For every x-slice and y-slice 
there is a recursively defined data structure. The reduction to rank space technique is applied on 
every 5 log log n-th level in the same way as in the Lemma [6l 

Given a query Q = [a, b] x [c, d] x (—00, e] we identify indices such that all cells C^, 

i\ < i < «2 and j% < j < 32 are entirely contained in Q. Then Q can be represented as a union of a 
query Q\ = [ciq, 60] x [cq, do] x (—00, e] and four (2, 1, l)-sided queries Q2 = [a, oq) x [c, d] x (—00, e], 
Q3 = (bo,b] x [c,d] x (-00, e], Q A = [ao,b ] x [c,c ) x (-00, e], and Q 5 = [a ,b ] x (d ,d] x (-00, e], 
where ao = x^, bo = Xj 2 _i, Co = y^, and do = yj 2 ~i- See Fig. [2] for an example. The query Q\ can 
be answered in constant time, and queries Qi, 1 < i < 5, can be answered using the corresponding 
x- and y-slices. Since queries Qi, 1 < i < 5, are equivalent to (2,l,l)-sided queries each of those 
queries can be answered in 0((loglogn) 3 + k) time. 

If the query Q is entirely contained in one x-slice or one y-slice, then Q is processed by a data 
structure for the corresponding x-slice resp. y-slice. Since the data structure consists of at most 
flog log n recursion levels, the query can be transferred to a data structure for an x- or y-slice at 
most v log log n times for v = 1^(2/(1+27)) ' H ence > tb- e total query time is 0(log log n + (log log n) 3 + 
k) = 0((loglogn) 3 + k). The space usage and construction time are estimated in the same way as 
in Lemma El □ 



Theorem 1 There exists a data structure that answers three-dimensional orthogonal range report- 
ing queries on the [U] x [U] x [U] grid in 0(loglogC7 + (log log n) 3 + k) time, uses 0(ralog 1+e n) 
space, and can be constructed in 0(n log 4 n log log n) time. 
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Proof : The result for the [n] x [n] x [n] grid directly follows from Lemma [7] and Lemma [TJ We can 
obtain the result for the [U] x[U]x [U] grid by applying the reduction to rank space technique |10|. 16] : 
We can use the van Emde Boas data structure [9] to find pred(e, S) and succ(e, S) for any e E [U] 
in O(loglogfJ) time, where S C [U] is P x , P y , or P z . Hence, the query time is increased by an 
additive term O (log log C/) and the space usage remains unchanged. □ 

Furthermore, we also obtain the result for d-dimensional range reporting, d > 3. 

Corollary 1 There exists a data structure that answers d- dimensional orthogonal range reporting 
queries in 0(log d_3 ra/(loglogra) d ~ 6 + k) time, uses 0{n \og d ~ 2+e n) space, and can be constructed 
in 0(n log rf+1+£ n) time. 

Proof: We can obtain a (i-dimensional data structure from a (d — 1) -dimensional data structure 
using range trees with node degree log £ n. See e.g. [2], [12] for details. □ 

Using Theorem Q] we can reduce the space usage and update time of the semi-dynamic data 
structure for three-dimensional range reporting queries. 

Corollary 2 There exists a data structure that uses 0(nlog 1+e n) space, and supports three- 
dimensional orthogonal range reporting queries in 0(logn(loglogn) 2 + k) time and insertions in 
0(log 5+£ n) time. 

Proof : We can obtain the semi-dynamic data structure from the static data structure using a variant 
of the logarithmic method [3 J . A detailed description can be found in [12] . The space usage remains 
the same, the query time increases by a O (log nj log log n) factor, and the amortized insertion time 
is 0(^^- log 1+£ n), where c(n) is the construction time of the static data structure. □ 

The result of Corollary [2] can be also extended to d > 3 dimensions using range trees. 

4 Three-Dimensional Emptiness Queries 

We can further reduce the space usage of the three-dimensional data structure if we allow 
O(loglogn) penalties for each point in the answer. Such a data structure can also be used to 
answer emptiness and one-reporting queries. As in the previous section, we design space efficient 
data structures for (2, 1, l)-sided and (2, 2, l)-sided queries. The proof is quite similar to the data 
structure of section [3] but some parameters must be chosen in a slightly different way. 

Theorem 2 There exists a data structure that answers three-dimensional orthogonal range re- 
porting queries on the [U] x [U] x [U] grid in 0(loglogC7 + (log log n) 3 + k log log n) time, uses 
0(nlogn(loglogn) 3 ) space, and can be constructed in time 0(ralog 4 n log log n) . 

For completeness, we provide the proof of Theorem [2] in the Appendix. Using the standard range 
trees and reduction to rank space techniques we can obtain a d-dimensional data structure for d > 3 

Corollary 3 There exists a data structure that answers d- dimensional orthogonal range reporting 
queries for d > 3 in 0(log rf ~ 3 n(loglogn) 3 + Hog log n) time, uses 0(n log d ~ 2 n(log log n) 3 ) space, 
and can be constructed in 0(n log d+1 n log log n) time. 
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Appendix. Proof of Theorem [2] 

Lemma 8 There exists a data structure that answers (2,l,l)-sided queries on the [n] x [n] x [n] 
grid in 0((loglogn) 3 + fcloglogra) time, uses 0(n(loglogn) 2 ) space, and can be constructed in 
0(n log 3 n log log n) time. 
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Proof: The data structure consists of the same components as the data structure of Lemma [6l But 
the size of x-slices and y-slices is reduced, so that each x-slice and each y-slice contains n 1//2 log p n 
points for a constant p > 2. The data structure Dt contains a point £ m m) f° r each cell 
Cij = XiHYj, Cij flP / 8, such that z m \ n is the minimal z-coordinate of a point in Cij fl P. 
The data structure Dt can contain up to n/log 2p n elements. Combining the results of Lemma Q] 
and Lemma [5j we can implement Dt in 0((n/ log 2p n) log n log log n) = 0(n) space, so that queries 
are supported in 0((loglogn) 2 + k) time. A list Ly contains all points in Cij sorted by their 
z-coordinates. For each x-slice X{, there are two data structures that support (1, 1, l)-sided queries 
open in +x and — x direction. For each y-slice Yj there is a data structure for (1, 1, l)-sided queries 
open in +y direction. For each x-slice and y-slice, there is a recursively defined data structure. 
As shown in Proposition 1 of [TT], the total number of elements in a data structure on the r-th 
recursion level can be estimated as s r (n) = 0{n l l 2r log p n^/log log n). The recursive sub-division 
stops when a data structure contains no more than log n elements. In this case, the data structure 
is implemented using e.g. the data structure of [2], so that queries are answered in O(loglogn) 
time and 0(logn(loglogn) 1+£ ) space. 

In the same way as in Lemma [61 the query Q can be represented as a union of (at most) one 
(2,l,l)-sided query on Dt, two (l,l,l)-sided queries on x-slices, and one (2,l,l)-sided query on a 
recursively defined data structure for a y-slice. Hence, the query time is 0((log log n) 3 ) if we ignore 
the time we need to output points in the answer. 

Unlike the data structure of Lemma [U we apply range reduction on every recursion level. Since 
the number of elements in a data structure on level r is s r {n) = 0{n l l 2r log p n\/\og log n), every 
element in a data structure on level r can be represented with log(s r (n)) = 0((l/2 r )logn + 
log log n) bits. Each data structure on level r uses 0(s r (n) log(s r (n)) log log(s r (n))) = 
0{s r {n) log(s r (n)) log log n) bits. The total number of elements in all data structures on level r is 
0{n2 r ). Hence, all level r data structures need 0(nlognloglog n + n2 r (log logn) 2 ) bits. Summing 
up by all recursion levels, the total space usage is Oin log n(log logn) 2 ) + ^^ x_1 n2 r (log logn) 2 
bits. The maximum recursion level r max = log logn + c r for a constant c r . Hence, the sec- 
ond term can be estimated as ^r=i x n ^ r 0°§ logn) 2 = 0(n log n(log logn) 2 ). If a data struc- 
ture on the recursion level r max contains m elements, then it uses 0(m(log log n) 1+e ) words of 
space because m < logn. All data structures on level r max contain O(nlogn) elements and use 
0(n log n(log log n) 1+e ) bits of space. Thus the data structure uses 0(n (log logn) 2 ) words of logn 
bits in total. 

The drawback of applying reduction to rank space on each recursion level is that we must pay 
a (higher than a constant) penalty for each point in the answer. Consider a data structure D r 
on the r-th level of recursion, and let P r be the set of points stored in D r . Coordinates of any 
point stored in D r belong to the rank space of P r . To obtain the point p G P that corresponds to 
a point p r £ P r we need 0(r) = O(loglogn) time. Hence, our data structure answers queries in 
0((log log n) 3 + k log log n) time. 

The construction time can be estimated with the formula 

c(n) = 0(n log 3 n log log n) + 2(n 1//2 / log p n)c(n x / 2 log p n) 

Therefore, c(n) = O (n log 3 n log logn). □ 

Lemma 9 There exists a data structure that answers (2, 2,1) -sided queries on the [n] x [n] x [n] 
grid in 0((loglogn) 3 + A; log logn) time, uses 0(n(loglogn) 3 ) space, and can be constructed in 
0(n log 3 n log logn) time. 
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Proof: The data structure is the same as in Lemma [8] but in each x-slice there are two data 
structures for (2, 1, l)-sided queries open in +x and — x directions, and in each y-slice there are two 
data structures for (2, 1, l)-sided queries open in +y and — y direction. 

The query is processed in the same way as in Lemma The space usage can be analyzed 
in the same way as in Lemma [8l Construction time can be estimated with the formula c(n) = 
0(n log 3 n log log n) + 2(n 1//2 / log p n)c(n 1 / 2 log p n) and c(n) = 0(n log 3 n log log n). □ 

Finally, we can apply Lemma [U and reduction to rank space and obtain the data structure for 
three-dimensional orthogonal range reporting queries on the [U] x [U] x [U] grid. 
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